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ABSTRACT 

M  This  paper  deals  with  the  structural  and  stability  properties  of  the 
averaging  approximation  scheme  for  linear  retarded  functional  differential 
equations.  Both  in  the  discrete-  and  in  the  continuous-time  case  the 
structure  of  the  approximating  systems  is  shown  to  be  analogous  to  the 
structure  of  the  underlying  retarded  equation.  Moreover,  it  is  shown  that  the 
approximating  systems  are  exponentially  stable  in  a  uniform  sense  if  the 
original  system  is  asymptotically  stable. 
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SIGNIFICANCE  AND  EXPLANATION 


Functional  differential  equations  play  an  important  role  in  a  large 
number  of  practical  problems  (wind  tunnel  control,  ship  stabilization,  steel 
mill,  chemical  engineering,  population  dynamics).  They  can  be  represented  as 
ordinary  differential  equations  in  infinite  dimensional  spaces.  Approximation 
methods  are  needed  in  order  to  solve  problems  like  dynamic  feedback 
stabilization,  optimal  control  or  parameter  identification.  For  some  of  these 
problems  Lyapunov  stability  properties  of  the  approximating  systems  are  of 
considerable  importance.  These  are  studied  in  the  present  paper  by  means  of 
an  analysis  of  some  fundamental  structural  properties  of  the  approximation 
scheme . 
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STRUCTURE  AND  STABILITY  OF  FINITE  DIMENSIONAL  APPROXIMATIONS 
FOR  FUNCTIONAL  DIFFERENTIAL  EQUATIONS 

Dietmar  Salamon 


1.  INTRODUCTION 

The  object  of  this  paper  is  to  present  some  new  results  on  the  averaging 
approximation  scheme  for  linear  retarded  functional  differential  equations  (RFDE) . 

The  averaging  approximation  scheme  has  been  invented  and  studied  by  several  Soviet 
authors  in  the  early  sixties  (see  e.g.  REPIN  [18)  i  further  references  and  a  detailed 
review  can  be  found  in  BANKS-BURNS  [21).  A  general  convergence  proof,  a  stability 
analysis  and  applications  to  optimal  control  problems  have  been  presented  for  the  first 
time  by  BANKS-BURNS  [1),(2].  Related  discrete-time  approximations  have  been  considered  by 
DELFOUR  [6),  REBER  [17],  ROSEN  [19).  Recently,  GIBSON  [9)  has  used  the  averaging  scheme 
for  approximating  the  solution  of  the  algebraic  Riccati  equation  associated  with  a 
retarded  system.  However,  there  remained  one  open  problem  in  the  convergence  proof  in  [9] 
which  has  not  yet  been  resolved.  This  is  the  question  whether  the  approximating  systems 
are  uniformly  exponentially  stable  for  sufficiently  large  N  if  the  underlying  RFDE  is 
stable.  In  [9]  this  has  been  stated  as  a  conjecture  without  proof.  We  show  in  section 
4.2  that  this  conjecture  is  in  fact  correct. 

Another  motivation  for  the  present  work  comes  from  some  recent  developments  in  the 
theory  of  retarded  systems  in  the  product  space  framework.  One  of  these  is  the 
introduction  of  so  called  structural  operators  for  the  state  space  description  of  RFDEs 
which  have  made  the  linear  theory  much  more  elegant  and  efficient  (see  e.g.  BERNIER- 
MANITIUS  (3),  MANITIUS  [141,  DELFOUR-MANITIUS  [7)).  They  have  led  to  a  number  of  new 
results  in  the  control  theory  of  RFDEs,  namely  on  problems  like  completeness  of 
eigenfunctions,  controllability,  observability,  and  the  linear  quadratic  optimal  control 
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problem.  Another  important  development  was  an  interpretation  of  the  adjoint  semigroup  in 


terms  of  the  underlying  REDE.  Extensions  to  neutral  systems  and  further  references  can  be 
found  in  SALAMON  [20  J. 

The  problem  has  not  yet  been  considered  whether  analogous  results  can  be  developed 
for  finite  dimensional  approximations  of  RFDEs,  in  particular  the  averaging  approximation 
scheme.  In  this  paper  we  fill  this  gap.  It  is  shown  that  the  approximating  systems 
satisfy  analogous  duality  relations  as  the  REDE  and  certain  structural  matrices  are 
introduced  which  play  an  analogous  role  for  the  approximating  systems  as  the  structural 
operators  do  for  the  REDE.  Moreover,  it  is  shown  that  these  matrices  actually  converge  to 
the  corresponding  structural  operators.  These  results  have  several  important 
consequences.  For  example,  they  lead  to  a  uniform  convergence  result  for  the  resolvent 
operators  and  they  are  crucial  for  the  proofs  of  the  stability  results  in  section  4.2. 

In  the  preliminary  section  2  we  give  a  brief  overview  over  some  recent  results  in  the 
theory  of  linear  retarded  systems  in  the  product  space  framework  and  describe  the 
averaging  approximation  scheme.  Section  3  is  devoted  to  the  study  of  the  structure  of  the 
approximating  systems  which  is  shown  to  be  analogous  to  the  structure  of  the  underlying 
REDE  under  several  aspects.  A  number  of  convergence  proofs  is  then  given  in  section  4.1 
and  two  stability  results  are  proved  in  section  4.2.  In  the  appendix  (section  6)  we  prove 
two  functional  analytic  results  which  are  frequently  needed  in  section  4.  In  particular, 

p 

we  give  a  quantitative  estimate  for  the  equivalence  of  I.  stability  and  exponential 
stability  for  strongly  continuous  semigroups. 

2.  LINEAR  RETARDED  SYSTEMS  AND  AVERAGING  APPROXIMATION 
2.  1  LINEAR  RETARDED  SYSTEMS 

We  consider  the  linear  retarded  functional  differential  equation 


x(  t )  ■  Lx^  ,  t  >  0, 


(2.1) 


-2- 


where  x(t)  e  R°  and  x^  is  defined  by  x^ix)  »  x(t+x),  -h  <  x  <  0,  h  >  0. 
Correspondingly  L  is  a  bounded  linear  functional  from  C "  C  [-h,0iRnl  into  Rn  given 
by 


0 

-  /  dn(x)$(T),  ♦  e  c, 

-h 


where  n(x)  is  an  nxn-matrix  valued  function  of  bounded  variation.  Without  lose  of 
generality  we  can  assume  that  n  is  normalized  which  means  that  n(x)  «  0  for 
x  >  0,  n(x)  »  n(-h)  for  x  <  -h,  and  n(x)  is  left  continuous  for  -h  <  x  <  0.  At  sosie 
places  we  will  assume  that  L  is  given  by 


L*  -  l  A  ♦(-h.)  +  /  A  ( x  )♦  ( T )  dx,  ♦  e  C, 
j=*0  -1  3  -h  U 


(2.2) 


where  0  -  hQ  <  . . .  <  hq  =  h  and  e  R  ,  j  «  0,  ...,  q,  as  well  as 
Aq  1  ( *  )  «  L2 (-h,0;Rnxn].  In  this  case  n  :  R  ♦  r"*"  is  clearly  given  by 


1(T>  -  -  V<— ,01<T)  •  V(-.-hj](T) 


-  J  AQ1(o)  do,  X  e  R. 


where  denotes  the  characteristic  function  of  the  interval  I. 

It  is  well  known  that  equation  (2.1)  admits  a  unique  solution 

x(  •  )  e  L2  t-h,«;Rnl  n  W.1'2  (0 ,«» ;I«n  ]  for  every  initial  condition  of  the  form 
loc  loc 


x(0) 


x!  x ) 


At), 


-h  <  x  <  0, 


(2.3) 


where  ♦  *  ( ♦ 1 >  €  R0  x  L2[-h,0iRn]  M2  .  This  solution  depends  continuously  on 
4  e  M2.  The  fundamental  solution  of  (2.1)  will  be  denoted  by  X(t)  e  R  ,  t  >  -h,  and 


-3- 


corresponds  to  the  initial  condition  X(0)  -  I,  X { t )  «  0,  -h  <  T  <  0.  It  can  also  be 
characterized  by  the  Volterra  integral  equation 

t 

X(t)  -  I  -  /  n(s~t)X(s)ds 
0 

and  its  Laplace  transform  is  given  by  MX)  *  where  A(X)  •  Xl-L(e*  ),  Jet,  is  the 
characteristic  matrix  of  (2.1). 

Proofs  of  these  facts  can  be  found  e.g.  in  HALE  [10]  or  DELFOUR-MANITIUS  [7). 

2.2  SEMIGROUPS  AND  STRUCTURAL  OPERATORS 

In  the  theory  of  RFDEs  -  as  well  as  of  other  types  of  integral  and  functional 

differential  equations  -  there  are  essentially  two  ways  of  introducing  the  state  of  the 

system  which  are  actually  dual  to  each  other.  The  state  of  system  (2.1)  in  the 

2 

"classical”  sense  is  the  pair  z(t)  »  (x(t),xt)  e  M  which  completely  describes  the  past 
history  of  the  solution.  Its  evolution  determines  the  strongly  continuous  semigroup 
S(t)  of  bounded  linear  operators  on  M2  defined  by 

s(t)$  -  (x(t),xt>  e  m2,  $  e  m2,  t  >  o, 

where  x(t),  t  >  -h,  is  the  unique  solution  of  (2.1)  and  (2.3).  The  infinitesimal 
generator  of  S(t)  is  given  by 

dam  a  -  {♦  e  m2^1  e  w1'2,  4>°  -  ♦ 1  (0 ) > , 

1  •  1 

A*  -  (L*  ,<P  ), 

where  w1,2  denotes  the  Sobolev  space  w1'2 [-h,0,Rn].  In  an  analogous  way  we  may 

2 

Introduce  the  semigroup  ST(t)  e  L(M  ),  t  >  0,  with  infinitesimal  generator  AT 
corresponding  to  the  transposed  RFDE 


-4- 


t  >  0 


(2.4) 


x(t)  -  L  xt  , 


Hie  duality  relation  between  (2.1)  and  (2.4)  can  be  described  by  means  of  an 

alternative  (dual)  state  concept  which  is  due  to  MILLER  (IS).  It  can  be  motivated  from 

the  fact  that  the  solution  of  the  RFDE  (2.1)  (t>0)  can  be  derived  from  the  initial 

1 

function  (t<0)  in  two  steps.  First  convert  the  initial  function  4  into  a  forcing 
term  of  suitable  length  which  determines  the  future  behaviour  of  the  solution.  Secondly 
determine  the  solution  which  corresponds  to  this  forcing  term.  The  dual  state  concept  is 
obtained  by  regarding  this  forcing  term  as  the  initial  state  of  the  system  rather  than  the 
solution  segment.  To  be  more  precise,  we  rewrite  equation  (2.1)  as 


c(  t )  =■  /  dn(t)*'t+T)  +  f- 1  (-t ) ,  x(0  )  -  f°. 


0  1  2 

where  the  pair  f  *  (f  )  e  M  is  given  by 


f°  »  4°,  f1 (a)  »  /  dn( T141 (t-o ),  -h  <  a  <  0. 
-h 


Now  the  initial  state  of  (2.5)  is  given  by  f  e  M  .  Correspondingly  the  state  at  time 
t  »  0  is  the  pair  w(t)  ■  (x!t),xt)  e  M2  where  x*"  e  L2[-h,0;Rn)  denotes  the  forcing 
term  of  the  shifted  equation  (2.5)  and  is  given  by 


xt(o)  -  /  dn(T)x(t+T-o>  +  f  (o-t),  -h  <  o  <  0. 


t  2  * 

The  evolution  of  this  state  (x(t),x  )  e  M  is  described  by  the  semigroup  ST<t)  (see 
e.g.  BERNIER-MANITIUS  [3]  or  SALAMON  (20)). 


Summarizing  this  situation,  we  have  to  deal  with  the  following  four  semigroups 


S  (t)  . 


-5- 


The  semigroups  on  the  left  correspond  to  the  RFDE  (2.1)  and  those  on  the  right  to  the 
transposed  RFDE  (2.4).  On  each  side  the  upper  semigroup  describes  the  respective  equation 
within  the  “classical"  state  concept  (solution  segments)  and  the  one  below  within  the  dual 
state  concept  (forcing  terms).  The  diagonal  relations  are  actually  given  by  functional 
analytic  duality. 

The  relation  between  the  two  state  concepts  can  be  described  by  means  of  so  called 

structural  operators.  These  have  been  introduced  by  BERNIER-MANITIUS  (3),  MANITIOS  (14), 

DELFOUR-MANITIUS  (7)  and  have  turned  out  to  be  a  very  elegant  and  efficient  concept  in  the 

2  2 

control  theory  of  RFDEs.  The  operator  F  8  l(n  )  maps  every  $  e  M  into  the 
corresponding  initial  state 


f*  ■  f  e  m 


of  (2.5)  which  is  given  by  (2.6).  The  operator  G  e  l(M  )  maps  every  forcing  term 
2 

fen  into  the  corresponding  solution  segment 


Gf  =  (x(h),x)  e  M 


of  (2.5)  at  time  h.  Thus  Gf  can  be  explicitly  described  as 


(Gf  ]  °  =■  (Gf)  ’(O), 

1  h  1 

(Gfl'lT)  «  X(h+T)f  *  J  X(h+T-s)f  (-s)ds,  -h  <  T  <  0. 
0 


Obviously,  G  is  bijective  as  an  operator  from  M  into  dom  A  and  its  inverse  is  given 


- 1  0  1 
(G  *1  -  *  (-h). 


(O)  *  J1  (-o-h)  -  /  dn(T)<fr  (T-0-h),  -h  <  o  <  0. 


i 


for  ♦  e  doin  A.  A  remarkable  fact  is  that  the  adjoint  operators  F  and  G  play  the 

same  role  for  the  transposed  equation  (2.4)  as  the  operators  F  and  G  do  for  the 
original  equation  (2.1).  Moreover,  the  following  result  has  been  proved  by  MANITIUS  (14) 
and  DELFOUK-MANITIUS  [7]. 

THE0RB1  2.  1 


(1) 

S(h)  -  GF,  ST(h)  *  FG. 

(ii) 

*  * 

FS(t)  *  ST(t)F,  S(t)G  *  GST(t),  t 

>  0. 

(iii) 

* 

If  ♦  e  dom  A,  then  F$  e  dom  A^ 

and 

A^F*  *  FA*. 

(iv) 

*  * 

K  fe  dom  At,  then  GA^f  »  AGf. 

We  close  this  section  with  a  concrete  representation  of  the  resolvent  operator.  For  this 
sake  we  introduce  for  any  X  e  t  the  operators  E^  :  (n  ♦  M2  and  T^:  M2  ♦  M2  by 
defining 


[E^x]0  “  x,  (Exx)  1  ( t )  *  e*Tx,  x  e  fn, 

(tx*)°  =  o,  [tx4>]  1  ( t )  =  /  eX(T"a  V<o)do,  *  e  m2. 

T 

Then  the  following  result  has  been  proved  in  MANITIUS  [13]  and  DELFOUR-MANITIUS  (7], 

PROPOSITION  2.2  Let  det  A ( A )  *  0.  Then 

(Xl-A)'1  =  ExA(X)_1Ey  F  ♦  Tx, 

•  -1  -1  *  * 

( Xl -A  )  =  FE.A(X)  Et  +  TT  . 

T  A  A  A 


l 


2.3  AVERAGING  APPROXIMATION 

In  this  section  we  briefly  describe  an  approximation  scheme  for  RFDEs  which  has  been 

studied  by  REPIN  [18],  BANKS-BURNS  [1],  [2],  GIBSON  [8]  and  many  others.  To  this  end  we 

N  2 

introduce  for  every  N  e  ■  the  linear  subspace  X  Cm  defined  by 


X1*  *  {  *  e  M2 1 <t> 1  (T 1  =  z.  e  Rn,  -  4  h  <  T  < 


•  .*0 


and  denote  the  corresponding  orthogonal  projection  by  p  :  M  ♦  X  .  This  subspace  can 

be  identified  with  Rn^N  ^  by  means  of  the  embedding  :  Rn^  ^  which 

T  T  T  n(N  +  1  ) 

associates  with  every  z  =  (Zg,...,z  )  8  R  the  pair 


.  N  ,  0 
[l  z]  =  z 


0' 


(.Nz]\t)  =  z.,  -  J  h  <  t  <  -  h,  j  -  1 . N. 


On  »n*N  *  we  always  consider  the  induced  inner  product 


,  s  T^N  „  „n(N+t> 

<  Z.w  >N  =  z  Q  w,  z.  w  e  R 


where 


X  0 


0  -  I 

•  N 


I 

N 


(2.8) 


The  corresponding  vector  and  matrix  norms  will  be  denoted  by  l»l  .  The  adjoint  operator 

N  ,  N  *  2  n  (N  +  1  ) 

»  «  ( 1  )  :  M  ♦  R  is  then  given  by 


[xN*)q  -  1"%]^  «  £  J  *1(T)dT,  j  =  1,...,N. 

-  -1  h 
N 

-8- 


I 


N  N 

Obviously#  the  operators  i  and  w  satisfy 


N  N  N  N  N 

w  l  -  id,  m  »  p  . 


(2.9) 


On 


#n(N+1 ) 


we  consider  the  differential  equation 


•N  N  N 

z  (t)  ■  A  z  (t)  ,  t  >  0# 


(2.  10) 


where 


(qW. 


N  N 
A0  A1 
r  -r 


(2.11) 


and 


-  lim  [TlfT  +  £)  -  n(T)l  ,  j  =  0,...,N. 

3  Tt-jh/N 


(2.12) 


N  N  -1  N  N 

In  an  analogous  way  we  define  the  matrix  At  *  (Q  )  Ht  where  the  A  ^  are  replaced  by 

NT  N  *  N 

( A  ^ )  for  j  *  0,1,..., N.  Then  the  adjoint  matrix  (A^,)  of  AT  with  respect  to  the 

inner  product  <  • , •  is  given  by 


<aV 

T 


(qn>-,<h y,  <„y  = 

T  T 


(2.13) 
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We  also  consider  the  differential  equation 


iN(t)  -  (A^)*wN(t)  ,  t  >  0, 
T 


(2.14) 


On 


Rn(N+1) 


The  following  theorem  has  been  proved  In  BANKS~BURNS  [2)  and  GIBSON  [9]. 


THEOREM  2. 3  Let  L  :  Q  *  Rn  be  given  by  (2.2).  Then  the  following  statements  hold. 
2  N 

( I)  For  every  $  e  M  we  have  4  =  lim  p  ♦  . 

(II)  There  exist  constants  M  >  1,  u  >  0,  such  that 


N 

\  < 


Me 


le  T  I 


<  Me 


for  every  t  >  0  and  every  N  e  *. 

2  2 

(111)  For  all  ♦  e  M  ,  f  e  M 

N  ,  N,  * 

„  A  t  .  „  (A  >  t 

s(t)4>  »  lim  t  e  *  ♦,  S  (t)f  *  lim  i  e  *  f 

N»“  N-»“ 

and  the  limits  are  uniform  on  every  compact  Interval  [0,T]. 

FULL  DISCRETIZATION 

A  fairly  general  and  extensive  study  of  full  dlscretitatlon  methods  for  RFDEs  can  be 
found  In  REBER  (17)  and  ROSEN  [19].  Since  the  aim  of  this  work  is  to  explore  the  special 
structure  of  the  averaging  approximation  scheme  described  above,  we  content  ourselves  with 
the  consideration  of  a  simple  one  step  Euler  approximation  for  the  ODE  (2.10)  which  has 
also  been  studied  by  DELFOUR  (6]  and  REBER  (17]  for  time  varying  systems 

Replacing  the  derivative  in  (2.10)  by  a  difference  quotient  with  step  sire  h/N,  we 


get  the  diffei  .e  equation 


(2. 15) 


N  h  N  N 

Vi  "  (I  +  N  A  >  V  “  > 


n(N+1 ) 

in  R  .  Since 


h  N  h  N 

i+nao  ?ai 


h  N 
N  AN 


(2.16) 


the  n(N+1 (-dimensional  1st  order  difference  equation  (2.15)  is  equivalent  to  the 
n-dimenslonal  (N+1)st  order  difference  equation 


N  .  N  N, 

h  'Vi  '  V 


N 

l 

j-0 


N 

Vi¬ 


le  >  0, 


(2.17) 


by  means  of  the  identification 


z 


N 

k 


VN  / 


e  r 


,n(N+1  ) 


(2.18) 


Equation  (2.17)  may  be  interpreted  as  a  direct  application  of  a  1-step  difference  approach 

to  the  RFDE  (2.1)  with  x^  approximating  x(khAI).  Finally,  note  that  this 

simplification  of  the  difference  equation  (2.15)  is  only  possible  because  of  the 

coincidence  of  the  step  size  h/N  for  the  time-discretization  with  the  mesh  size  of  the 

N  2 

spatial  discretization  in  the  subspace  X  c  M  . 


3.  THE  STRUCTURE  OF  THE  APPROXIMATING  SYSTEMS 

It  is  the  goal  of  this  section  to  analyse  in  detail  the  structure  of  the 
approximating  systems  (2.10),  (2.14)  and  (2.15)  respectively  (2.17).  It  is  shown  that 
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I 


there  is  a  strict  analogy  to  the  structure  of  the  underlying  RFDE  (2.1)  as  it  has  been 


M 

described  in  section  2.2.  In  particular,  there  are  certain  structural  matrices  F  and 

G1*  playing  the  same  role  for  the  approximating  systems  as  the  operators  F  and  G  do 
for  the  RFDE  (2.1). 


3.1  THE  STRUCTURAL  MATRICES 

Starting  from  (2.17),  we  observe  that  there  is  another  way  of  transforming  this 
(N+1)st  order  difference  equation  into  an  equivalent  1st  order  equation.  For  this  sake 
let  us  rewrite  (2.17)  as 


N  ,  N 

h  (*k 


+  1 


N. 

V 


k 

l 

j-0 


N 

V 


+  f, 


k+1 


k  >  0, 


(3.1) 


N  N 

where  A^  :•  0,  f^  s-  0  for  j  >  N  and 


N  N 

f.  ■  X.  , 
0  0 


.N  V  .»  N  , 

'« ’  ,L  w 1  ■  '• 


..,N 


The  forcing  term 

i  ^ 

fN  -  : 

\  <} 


e  r 


n(N+1 ) 


(3.2) 


may  be  considered  as  the  initial  state  of  (3.1)  since  it  contains  all  the  information 
which  is  needed  for  determining  the  future  behaviour  of  its  solution  ,  k  >  0. 
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where 


Correspondingly  the  state  at  instant  k  e  H  is  given  by  wj|  e  Rn*N  ^ 


N  N 

Vo  '  V 


N  ’'+r"1  N  N  N 

V*  Aj  Vt-l-j  fk+i'  1 


(3.3) 


.,N 


Then  it  is  easy  to  see  that  w^  satisfies  the  1st  order  difference  equation 


V,  ■  (I  +N  V  k  >  °‘ 


(3.4) 


since 


1  +  £  <4>‘ 

N  T 


h  N 

1 4 


S- 


Note  that  (3.4)  can  be  regarded  as  a  one  step  Euler  approximation  for  the  ODE  (2.14). 

We  conclude  that  there  are  two  state  concepts  for  the  difference  equation  (2.17), 
namely  (2. IB)  and  (3.2-3),  both  of  which  lead  to  a  first  order  difference  equation  in 
r"(N+1>,  namely  to  (2.15)  and  (3.4).  The  relation  between  these  two  state  concepts  can 
be  described  by  certain  structural  matrices  F^1  and  G?1.  Before  defining  these  matrices, 
we  introduce  the  concept  of  a  fundamental  solution  for  equation  (3.1). 


DEFINITION  3. 1  The  fundamental  matrix  of  equation  (3.1)  is  the  sequence 
XN  e  Rnxn^  R  >  0>  deflned  by 


S  (XN  -  XN) 
h  1  k+1  k' 


„  N  N 


3-o 


Vi'  k  e  ■'  xo 


(3.5) 
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Remark  3.  2 


(1)  By  induction,  it  is  easy  to  see  that 


f*N  -  vN,  .  y  .N 

Xk+1  V  ^  Xk-j  Aj 


N  ,„N 
h 


,  k  e  n. 


(3.6) 


(ii)  The  solution  of  (3.1)  is  given  by 


k-1 

N  h  r  N  ,N  ,  „  „ 


(3.7) 


Now  we  introduce  the  matrices 


I  0 
°  *N 


0  A" 


(3.8) 


and 


N-N,  ^ 


K  0 


XN  '  ' 


(3.9) 


N 


Then  it  is  easy  to  see  that  f 


N  N 

F  z0  If 


C  Rn(N  1)  ^  defined  by  (2.18)  and 


fN  e  Rn*N  1  *  is  the  forcing  term  of  (3.1)  defined  by  (3.2).  Moreover,  if  xj^,  k  >  0, 

N  n(N  +  1) 

zN  e  R  is  defined  by  (2.18),  then  it  follows  from 


is  the  solution  of  (3.1)  and 


remark  3.2  that  z^  =  GNfN.  Making  use  of  these  facts,  one  can  easily  establish  the 
following  result  which  is  strictly  analogous  to  theorem  2. 1.  The  proof  is  left  to  the 
reader . 
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PROPOSITION  3.3 

(1)  (X  ♦  |  aV  -  gV,  (I  ♦  |  (a")V  -  rV. 


(U)  fV  -  <a")V  .  angn  -  gN(a")\ 

(1,n  NANt  (AT>\n  aVn  « 

( 111 )  Fe  *  e  F  ,  e  G  *  G  e  #t>0. 


(lv) 


(gV1 


M 

O 

• 

• 

• 

o 

1 _ 

o 

• 

• 

o 

_J 

•  • 

•  •  N  _  N  _ 

•  •  *  N 

•  •  —  T  «  —  T 

0  •  h  h 

• 

•  •  •  0 

•  • 

•  •  0 

•  •  • 

•  •  •  • 

0  *  »N  .  .  . 

A0  ^*-1 

_  — 

L_  —* 

N  N  N  N 

We  conclude  that,  for  any  solution  z  (t)  of  (2.10),  the  function  w  (t)  ■  F  z  (t) 


satisfies  (2.14)  and,  conversely,  for  any  solution  wN(t)  of  (2.14),  the  function 
rN(t>  -  GNwN(t)  satisfies  (2.10). 


3.2  SPECTRAL  THEORY 

In  this  section  we  give  a  brief  overview  over  some  spectral  properties  of  A11  and 

N  * 

<At)  which  are  analogous  to  wellknown  results  in  the  theory  of  RFDEs.  In  particular, 
we  will  see  that  the  rational  complex  nxn-matrlx  valued  function 


Ax,  -  Xl  -  LN(X),  ln,x, 


N 

I 

j-0 


AN  (— - )’ 

"  j  N+Xh 1 


(3.10) 


plays  precisely  the  same  role  for  the  approximating  systems  as  the  characteristic  matrix 
MX)  does  for  the  underlying  RFDE  (2.1).  Moreover,  we  introduce  the  matrices 


E 


N 

X 


I 

N 

N+\h 


<-JJ-)N  i 
lN  +  Xh' 


e  t 


n(N  +  1  )xn 


(3.11) 
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I 


and 


rN  =  & 

X  N 


N+Xh 


I  • 


•  • 

•  • 

•  • 


0  W  1  •  • 


0 
N 

N+Xh 


e  t 


,n(N+1  )xn(N+1 ) 


(3. 12) 


for  X  e  X  f  -  N/h. 


LEMMA  3.4  Let  )  f  t,  X  f  -N/h,  and  z,w  e  (fn^N  ^  be  given.  Then  (Xl-A*Sz  *  w  If  and 
only  If 


N  N 

2  =  EX20  +  V 


(3. 13) 


and 


.N , ,  ,  .  N  T  N  N 

A  (X  )*Q  »  (E^ )  Q  F  w  . 


(3. 14) 


N  N  N  N 

PROOF  Clearly  (Xl  -  A  )z  =  w  if  and  only  if  (XQ  -  H  )z  “  p  w  or  equivalently 


Xz0 '  E  Vj 

j-0  J  ] 


(3.15) 


z.  =  rrrr  (::  w  +  z  i ,  j  -  i, . . . ,n. 


j  N+Xh  N  j  j-1 


(3.16) 


Equation  (3.16)  is  equivalent  to 


(_N_,jz  +  h  y  ( — S-, 

'm  +  XVi '  zn  «  i  'u+lh' 


j  N  +  Xh  0  N  +  ,  N+Xh  j+1-v 

■>  v»  1 


3  “  1,  .  .  .  ,N 
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*nd  hence  to  (3.13).  If  this  is  satisfied,  then  (3.15)  is  equivalent  to 


AH(X)z„ 


N  r  N  r  h  ?  , _ N  ,  V 

*  *0*0  -  ^  Aj  N+Xh  ‘ 


j+1-vJ 


,  h  p  ,  N  V  p  N 
W0  N  N+Xh’  j  “j+1-V 


,  N  TNN 
(E^)  0  F  w. 


Note  that  the  above  lemma  is  strictly  analogous  to  a  well  known  result  in  the  theory  of 
RFDEs  (see  e.g.  HALE  (10),  DELFOUR-MANITIUS  I?]).  It  has  several  important  consequences 
which  are  summarized  in  the  proposition  below  and  can  be  proved  straight  forward. 
Statement  (i)  can  be  found  in  BANKS-BURNS  [2).  Statement  (iv)  is  the  analogon  of 
proposition  2.2. 


PROPOSITION  3.5 

(1)  Let  1  e  t,  X  f  -N/h,  then  X  e  (J(AN)  if  and  only  if  det  AN ( X )  ■  0. 

( ii)  X  -  -N/h  e  ct(AN)  if  and  only  if  det  A^j  •  0. 

(iii)  «((»;>*)  -  o(a”)  -  c(aN). 

(iv)  If  X  (■  -N/h  and  det  A  (X)  ^  0,  then 


(XI-aV1  .^(xr'iE^Vp8*^ 


(XH^V’-Ayaf'^iVt  (r”)T. 


REMARK  3.6  A  solution  x(t)  of  the  RFDE  (2.1)  is  said  to  be  small  if  it  vanishes  after 
some  finite  time  T  (HENRY  [11)).  If  L  t  *  «"  li  given  by  (2.2)  and  if  AQ 1 < • )  1  0, 
then  there  exist  nonzero  small  solutions  of  (2.1)  if  and  only  if  det  A^  «  0  (MANITIUS 
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[14].  Now  not*  that  for  sufficiently  large  N  this  means  that  dst  A^  ■  0  and  hsncs 
-N/h  e  o(AN)  (proposition  3.5  (11)).  This  indicates  that  ths  generalized  elgenmodes  of 
(2.10)  rsspectivsly  (2.14)  corresponding  to  ths  eigsnvalus  X  *  -N/h  play  the  role  of 

the  small  solutions  in  the  approximating  systems.  Moreover,  note  that  the  solutions  of 

U 

the  difference  equation  (2.15)  starting  with  generalized  eigenvectors  of  A 
corresponding  to  X  *  -N/h  are  precisely  those  solutions  which  vanish  after  a  finite  time, 

4.  CONVERGENCE  AND  STABILITY 

Having  introduced  a  number  of  operators  for  the  approximating  systems  which  are 
analogous  to  well  known  operators  in  the  theory  of  RFDEs,  we  may  pose  the  question,  if  - 
and  in  what  sense  -  these  operators  converge,  ttiis  problem  will  be  considered  in  the  next 
section. 

4. 1  CONVERGENCE 

We  begin  with  some  preliminary  facts. 


REMARK  4. 1 

(i)  It  is  easy  to  see  that  the  function  n**  :  R  ♦  R1'*1'  defined  by 


N  V-  1  V  V. 

n ( T )  -  Urn  n  ( T ) ,  —  h  <  t  <  —  ,  k  e  Z, 
xtkh/N  N  N 


satisfies  the  inequality 


/  I  n  ( r )  -  nN  (  t  )  I  dr  <  jj  var(h)  . 
-h 


(4.1) 


( 1 i )  For  every  k  e  t,  k  f  -N/h,  let  us  define  the  function  e.  :  [— h , 0 )  ♦  t  by 


ex(T)  *  -  n  h  4  T  <  -  Vh'  j  ■  0 . N- 
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Then  it  la  well  known  that  the  limit 


lim  sup  |e*T  -  e*?(T)|  *  0  (4.2) 

N-m*  -h<t<0 


exists  uniformly  on  bounded  subsets  of  the  complex  plane. 


The  following  convergence  result  for  AN(X)  has  been  shown  by  BANKS  and  BURNS  [2]. 
For  completeness,  we  present  an  alternative  and  simplified  proof. 


LEMMA  4.2 

(i)  A  (X)  converges  to  A(X)  uniformly  on  every  bounded  subset  of  the  complex  plane. 

|| 

(11)  For  every  a  >  0  there  exists  a  constant  cq  >  0  such  that  |L  (X)|  <  cq  for  every 
NeN  with  N  >  ah  and  every  X  e  $  with  Re  X  >  -a. 

PROOF  Note  that 

L(e*  )  *  /  e*Tdn(r)  *  -n(-h)e  -  X  /  n(T)e^TdT 

-h  -h 


and,  by  (3.10)  and  (2.  12), 


LN(X)  -  l  [nV  h)  -  nV  i  h)]  (^)j 


j»0 


N  N  N  .hr  N  ±li  .  ._N_.j 

n  (-h)  W,’  X  N  n  (  N  M  N-tXh 


0 

n  <-h)  e^(-h)  -  X  /  nN  (T )  ej(t)  dr  . 
*  -h 


Thus  statement  (i)  follows  immediately  from  remark  4.1.  Statement  (li)  follows  from 

N  N 

(3.10)  with  c  -  VAR(n)  sup  {<~ — r>  IN  >  ah)  <  “  . 

<1  N-an 
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For  the  next  result  we  need  the  space 


M  -  Rn  x  L  (-h,0($n) 

0  1  • 
endowed  with  the  norm  l$l  -  max  (!♦  I ,  !♦  I  }  for  $  e  M  . 

M*  L* 

THEOROI  4.3  The  limits 

lim  IE,  -  ,Vl  “  0  ■  lim  IT,  -  ,VV. 

N*<*  A  N-»  HM2 ,h“ ) 

exist  uniformly  on  bounded  subsets  of  the  complex  plane. 

PROOF  The  statement  on  E^  Is  an  immediate  consequence  of  remark  4.1  (ii),  since 
iNE^x  *  (x,e^(*)x)  e  M  for  x  e  $n,  N  e  ■  and  \  e  X  t  -N/h. 

In  order  to  prove  the  second  part  of  the  theorem,  let  us  first  define  e(X,r)  e  M 
by 

e(X,T)°  -  0,  e( X,r) 1 (o)  - 

for  X  e  <t  and  -h  <  T  <  0.  Then 


0  ,  -h  <  a  <  T, 

-X(T'0)  ,  t  <  o  <  0, 


[TxpN*l  1  (T  ) 


<  e(X,T),  p  ♦  >  =  <  p  e(X,T),$  > 


2 

for  all  *eM,Xet,Ne«  and  T  e  (-h,0J.  Moreover,  the  closure  of  the  set 

2  N 

(e(X,T)  |  | X |  <  c,  -h  <  T  <  0}  in  M  is  compact  and  thus  p  converges  uniformly  on 

N  2  «• 

this  set.  Hence  T^p  converges  to  T^  in  i(M  ,M  )  uniformly  on  bounded  subsets  of 
the  complex  plane . 

Secondly,  note  that 


,  N  N  , 
[l  T.zl 


(T) 


j 

h  r  _ N_.v 

N  L ,  N+Xh'  Zjt1-V 
v-1  J 


/  e^(o)  [  iNz]  1  (-  -a)  do 

-jh/N 
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I 


and  hence 


,  N.,  1  .  .  N_  N  N, ,  1 ,  .  . 

|  (T^p  *]  (X)  -  lI  Tx»t  4)  ( T )  | 


<  I  [T.p%] 1  (T  )  -  fr-pN*)1(-  ^)|  +  /  |e*°  -  Oj  (O)  |  |  [pN41  *  (-  jp1  -  o)  |  do 

A  A  N  -jh/H  *  N 

for  -  ^  h  <  T  <  -  h,  j  «  1,...,N.  Thus  the  statement  of  the  theorem  follows  from 

N  N 

remark  4.1  (11)  together  with  the  fact  that  the  set 


{ rrx41 1 1  |X |  <  C,  ♦  e  m2,  141  <  1} 


Is  equlcontlnuous  for  c  < 


THBORP1  4.4  For  every  4  e  M 


N  N  N, 

F4  “  lim  l  F  w  4  , 


PROOF  We  prove  this  result  In  3  steps.  The  first  step  is  a  formula  for  the  operators 

N  ^  N  N 

11  F  and  F  *  . 


Step  1  (»NF4]q  «  [FN»N410  ”  4°  and  for  j  * 


[wNF41  "  *  r  J  n(x-j4i)  [♦1(t)  -  ♦’(T-h/W)]  dx 
3  -h 

(F^s^4l  j  “  -  jj  /  nN(x-^h)  [♦’(x)  -  *1(x-hAl)]  dx 


Proof  Let  us  define  4  (x)  :  -  0  for  x  4  (-h,0).  Then 

-JTh  0 

[»KF41  ”  f  /  J  dnfx^’lx-c)  do 

Jh  'h 

N 


I 


do  -  ^  /  n(T)  [*1(T+^tl)  -  ♦1(t+'^-4i)] 

-h 

-  ♦  ’ (T-jJ> ]  dT 

h/N 

-  h  J  n(T-^h)  [♦  <  T )  -  ♦’<T-|>]  dT 

-h4h 


r  I  dn ( t )  /  *  <c 


-h 


-  ~  n(-h)  /  4> 1  < a > 

N 

-h*jh 

-  ^  n(-h>  /  [♦ 1 ( t ) 


h/N 


^  J  n(T-jJh)  [♦’  <  T  ) -4> 1  (T-^)  ]  dT 
“h 


and 


IfV^I  -  [  A  [,%] 

j  L  .  v  v-  i 

J  v»3  J 


j*1 


N  N 


I  !  [hN<T-^pi)  -  nN  (  T-jjh )  ]  ♦ 1  ( T )  dT 


v-  j  v—  i+ 1 


^  i  [hN(T--l^h)  -  hN(T-^h)]  ♦ 1  <  t  )  dT 
~h 


=  -  ~  j  hN(T-^h ) [♦ ' ( T )  -  ♦’(T-^)]  dT 


-h 


for  j 


1 ,  •  .  • , N •  Tliis  proves  step  1. 


step  2  ,fN,n  <  SUP  O#  VAR  { n ) }  V  N  e  *. 

Proof  By  the  well  known  convolution  inequality,  we  have 


dT 
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N  „  ,  N  N 

I  I  [FNz]  ,|2  *  l  I  [  A 


J-1 


j-1  v-j 


u  v-j+1 


N 

l 

v-1 


N 


<  [  i  i\ir  i  i*j- 

j-i 


<  [  var ( n ) ] 2  l  I*  r 

j-1  5 


and  hence 


,  N  ,2  .  .  2  h  r  .  .  N  ,  .  2 

,F  Z,N  *  Z0  +  H  £,  ltF  Zlj' 


<  U0|2  ♦  [VAR(n>]2  t  l  |Z  f 
j-1 


<  [max  {1,  VAR(n)}]2lzi2 


for  z  e  Rn^N  ^  and  New.  This  proves  step  2. 


N  N  N  2 

Step  3  F$  -  Lim  I  F  M  V  $  e  M  . 

N  ♦<* 

Proof  Let  us  first  assume  that  is  continuously  differentiable  and  that 

^ 1 ( 0 )  -  0.  Moreover  let  us  define  $  ^ ( T  >  -  0  for  T  >  0,  Then  it  follows  from 

step  1  and  (4.1)  that 


,  ,  N  ,  ,  .  N  fL  ,  . 

I  [»  F*]  -  [F  »  ♦!  | 


h/N 

/  [n  <t-^>  - 

-  ♦1(Th/S>]  di| 

-h 

■°  N 

In  (t)  -  n(r)|  dx 

N  , 

sup  .  1 

/  i  ( o )  do  | 

•h 

■h<r<h/N 

i-h/N 

h  *1 

<  -  Var ( n )  I*  I 

N  m 
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and  hence 


.  „ ,  NN  N 
IF*  -  l  F  »  $1 


<  IF*  -  pNF*l  +  lir^F*  -  A%l 


<  IF*  -  pNF*l  +  M  1  |(hNF*]  -  [FNnN*]  |211/2 

j«1  :  1 


<  If*  -  p  f*I  +  ^ 


3/2 


VAR(n)  i*1! 


For  any  ♦CM  the  statement  follows  from  the  Banach-Steinhaus  theorem  and  step  2.  This 
proves  the  theorem.  ^ 


Combining  the  above  convergence  results  with  the  concrete  representation  of  the  resolvent 
operators  given  in  proposition  2.2  and  proposition  3.5,  we  obtain  the  following  result. 
The  proof  is  a  straight  forward  application  of  lemma  6.1  and  will  be  omitted. 

COROLLARY  4.5  The  limits 

lim  l(Xl-A)'1  -  lN(Xl-AN)“’itNl  2  -  0  , 

N~»  L(M  ,  m”) 

*  - 1  n  N*-1N 

lira  I(XI-At)  -  1  (Xl-(AT)  )  it  I  -  0 

Nx»  L(M^) 

exist  uniformly  on  those  bounded  subsets  of  the  complex  plane  which  are  uniformly  bounded 
away  from  the  zeros  of  det  A ( X ) . 

THEOREM  4.6 

lim  IG  -  TNGNxtl|  2  «  -  0  * 

Nx»  1(M  ,M  ) 


4- 


PROOF  We  establish  this  result  In  three  steps- 

Step  1  Let  Xy  j  >  0,  be  given  by  (3.5)  and  let  us  define 


XN(t)  *-X;.jh<t<^h.J.  0,1.2 . 


Then  X**(t)  converges  to  X(t)  uniformly  on  every  compact  interval  [0,T] 
Proof  For  every  ken 


/•nT  tx"  -  x"l 

k  V+1  v 

v-0 


k- 1  v 

T  +  tL  V  V  aN  *N 

1  N  i  Av-j  *j 

v«0  j»0  J 

„  X"1  k_1  U  M 

i  + ^  y  an  xn 

N  *  L  v-1  1 

1*0  v-1  J  J 


n  a  n  ‘nr™  xj 

j»0  J 


T  N,  iiii  . 

*  I  -  ]  n  ( a-~)  x  ( a)  de 

o  * 


k  k-M 

and  hence  for  ~  h  <  t  <  — r-  h 

N  N 


X(t>  -  XN(t)  -  X(t)  -  X(^) 


kh/N 


♦  /  [nN(s-~k>  -  his-—1)]  X(s)<5b 


kh/N 

*■  /  n  (s-^lfx  (s)  -  X (  • )  ] da ■ 

0 


Thus  the  desired  convergence  result  follows  from  Gronwall's  lemma. 


Step  2  Let  z  e  Rn(N  +  1  '  and  j  e  {! . N) .  T^en 


[>NGNz]^  “  XN(h)z  +  /  XN  (s)  [  lNz]  1  (a-h)ds, 
0 


[iNGNzl  (t)  »  XN(h+r)z  +  /  XN(s)[vNz]'(  s-h+^)da,  -  ^  h  <  T  <  -  h. 

u  Q  N  N  N 


Proof  If  -  h  <  T  <  -  h,  then 
N  N 


N-j-1 


1+1. 


XN(h+t)z  +  £  /  XN(s)  [  iNz]  1  ( a-h+^1  da 


j. 


XN(hH)z  +  /  XN(a)  [iNz]  1  (s-h+jj^)  da. 
0  o 


In  the  case  j  *  0  this  equation  leads  to  the  desired  expression  for  [i^G^x]^. 


Step  3  lim  IG  -  lNGN*NI  „  *  0* 

N~  L(M2,m“) 

Proof  First  note  that  the  functions  [Gf]  e  C,  f  e  M  ,  Ifl  <  1,  are  equi continuous 

1  2 

aince  the  canonical  embedding  of  W  '  into  C  is  a  compact  operator. 

,  n(N  +  l ) 

Now  let  z  e  R  .  Then,  by  step  2, 


N  N  N  ,  0 
(Gl  z  -  l  G  z] 


[x(h)-XN(h)  ]z  +  /  [x( a )  -  XN(s>)  t iNz! 1 ( s-h)  da 


and  for  -  jj  h  <  T  <  -  h,  j  =  1,  ....  N, 
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This  has  also  been 


and  moreover  that  this  convergence  is  uniform  for  bounded  f  e  M  . 

proved  by  REBER  (17,  theorem  7.5]  under  the  condition  that  L  :  C  ♦  Rn  is  given  by  (2.2). 

N  2  N*  2 

Let  us  now  introduce  the  operator  families  S  (t)  e  i-(M  ),  S  (t)  e  HM  ),  t  >  0,  by 


sN(t) 


Nf  h  N 

‘  t1  +N  A 


sr«t> 

T 


-  ."[* 


h  ,  N  ‘,k  N 

-  <Am>  J  *  , 

n  T  ‘ 


|  <  t  c  ^  h,  k  -  0,1, 2,. 


(4.3) 


Then  the  following  result  is  a  direct  consequence  of  theorem  4.4  and  theorem  4.6  together 
with  the  factorization  results  (theorem  2.1  (i)  and  proposition  3.3  (i)). 

COROLLARY  4.7 

(I)  For  all  *  e  M2,  fen2 

S(t>*  -  lim  SN(t)$,  S*(t)f  *  lim  S^*<T>f 
Nx»  N»” 

and  the  convergence  la  uniform  on  every  compact  Interval  [0,T). 

( II)  For  every  k  e  N 

11m  IS(kh)  -  SN(kh)l  ^  -  11.  I  S*(kh)  -  S^*(kh)l  ,  -  °. 

n-  1(M  ,m”)  nx»  Km  ) 

PROOF  It  only  remains  to  note  -  for  the  proof  of  statement  (11)  -  that,  by  lemma  6.1, 

lim  I  S ( h )  -  GiNFNaNl  2 
Nx»  l(«  ,  M  ) 

•  *  N  N  T  N  * 

-  lim  IF  C.  -  l  (F  )  w  G  I  «  *  j  *  O'  ■ 

NX”  i(  (M  )  ) 
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Stateaant  (ii)  of  the  above  result  la  apparently  new.  The  strong  convergence  of  statement 
( 1)  has  bean  stated  without  proof  by  DELFOUR  [6].  The  strong  convergence  of  SN(t)  has 
been  shown  by  REBER  [17]  and  ROSEN  [19]. 

4.2  UNIFORM  STABILITY 

It  Is  a  simple  consequence  of  corollary  4.7  that  the  discrete  time  systems  (2.15)  and 
(3.4)  are  stable  in  a  uniform  sense  if  the  underlying  RFDE  (2.1)  is  exponentially 
stable.  More  precisely,  we  have  the  following  result. 


THEORPI  4.8  Let  w  <  0  and  suppose  that  det  A(X)  +  0  for  every  let  with 
Re  X  >  w.  Then  there  exists  an  e  N  and  a  constant  Y  >  0  such  that 


I 


<  Ye 


ukh/N 


for  every  N  >  N^. 

PROOF  It  follows  from  a  well  known  result  in  semigroup  theory  that  there  exists  a 

wk0h 

k  e  *  euch  that  IS(k  h ) I  <  e  .  By  corollary  4.7  (li),  this  implies  the 

L(H  ) 

existence  of  an  NQ  e  H  such  hat 


.  „  k  N  ink  N 

KI  +  “•  *N)  °  *N  <  e  ,  N  >  NQ  . 


Moreover,  it  follows  from  corollary  4.7  <i),  that 


rr  y,  m  # 

Y  e  sup  (KI  +  $  A  >  «  |l  -  0,  ...,  knN-1,  N  e  ■}  <  -. 

N  N  U 


We  conclude  that  the  following  inequality  holds  for  N  >  and  k  “  vk^N  +  *  with 
v  e  m  and  t  e  {0,  . ..,  kQN-1 } 
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I 


1(1  +  £  AN)kiN  <  »«  ♦  £  aVinki  +  jj  an)  0  |V 


uikjjh  (DVkgh 

<  Ye  e 


<  Ye 


<o(vk0+t/N)h 


re 


tokh/N 


It  follows  easily  from  lemma  4.2  that  the  stability  of  the  RFDE  (2.1)  also  implies 
the  stability  of  the  approximating  continuous-time  systems  (2.10)  and  (2.14)  if  N  is 
sufficiently  large  (the  precise  arguments  are  given  in  the  proof  of  theorem  4.9  below). 
However,  a  uniform  estimate  in  the  spirit  of  theorem  4.8  has  not  yet  been  proved  in  the 
literature  on  these  approximation  schemes.  It  has  been  stated  as  a  conjecture  by  GIBSON 
[9]  and  provides  -  in  that  paper  -  a  crucial  step  in  the  convergence  proof  for  the 
solutions  of  the  algebraic  Riccati  equation.  REPIN  (18)  also  claims  the  uniform  stability 
of  the  approximating  systems  (2.10),  however,  his  arguments  are  extremely  unclear  and  it 
seems  almost  impossible  to  convert  them  into  a  rigorous  proof.  The  following  theorem 
closes  this  important  gap  in  the  approximation  theory  of  RFDEs  and  may  be  considered  as 
the  main  result  of  this  paper. 


THEOREM  4.9  Let  L  :  C  ♦  h"  be  given  by  (2.2)  and  let  the  RFDE  (2.1)  be  exponentially 
stable .  Then  the  approximating  systems  (2.10)  and  (2.14)  are  uniformly  exponentially 
stable  for  sufficiently  large  N.  This  means  that  there  exists  an  NQ  e  *  and  constants 
£  >  0,  y  >  1  such  that 


N 

•  A  t, 
le  I 


N' 


(A\ 

•e  T  I 


<  Ye 


for  every  t  >  0  and  every  N  >  . 
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N  *  u 

PROOF  First  not*  that  the  statement  on  (A^,)  follows  from  that  on  A  .  Secondly,  It 
follows  from  theorem  6.2  and  the  exponential  estimates  in  theorem  2.3  (ii)  that  it  is 
enough  to  show  that  there  exists  an  e  H  and  a  constant  c  >  0  such  that 

«  N 

/  le*  ‘zl2  dt  <  c2lzl2  (4.3) 

0 

for  every  z  «  Kn(N+1'  and  eVery  N  >  Nfl.  We  will  prove  this  in  5  steps. 

Step  1  There  exists  an  nq  e  B  such  that  det  A  (1)  f  0  for  every  1  e  (  with 
Re  X  >  0  and  every  N  >  . 

jq 

Proof  By  lemma  4. 2 ( ii ) ,  the  complex  function  det  A  (1)  cannot  have  a  zero  in  the  closed 
right  halfplane  outside  the  disc  of  radius  VAR(n)  centered  at  the  origin.  Inside  this 
disc  the  nonexistence  of  unstable  eigenvalues  of  A  follows  frost  lenxaa  4.2(i)  if  N  is 
sufficiently  large. 


Step  2  For  New  let  us  introduce  the  matrix 


-1 


e  r 


_NxN 


Then  there  exist  constants  cq  >  ' 


'•'"‘'I*.  *  V  ^  V  *  »  °  V  "  e  ■* 


such  that 


T  N  N 

Proof  rirst  of  ail  it  la  easy  to  see  that  x  a  x  *  0  for  every  x  e  R  and  every 
N  €  ■.  Hence  it  follows  front  a  well  known  result  in  semigroup  theory  that 
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v  t  >  o  v  n  e  * 


I^NxN  *  1 


Moreover 


a\ 


N 

h  1 


Nt/h 
1 1 


(Nt/h) 


N-1 


•  • 


(N-1) ! 


Wt/h 
1 1 


and  hence 


N  N~  1 

at,  .  r  1  ,Nt.k  -Nt/h 

B  NxN  «  Jn  k!  '  h>  6 
k«0 


for  every  t  >  0  and  every  N  e  ■.  Since 


k  -Nt/h  ..  .  ,  ^.k+1 

t  e  dt  »  kl  (h/N)  , 


this  implies 


[  '^'nxN  dt  4  Ja  5  "  h 
o  k=o 


N-1 


Together  with  (4.4)  this  estimate  proves  the  statement  of  step 
More  precisely,  €g  >  0  may  be  chosen  to  be  any  constant  less 


n(N+1 ) 

Step  3  For  every  z  e  R  and  every  N  e  ■ 


/  It"  zl*  dm  <  — -  Izl*. 
'  ^  tU)  N  tg  N 
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2  (theorem  6.2,  p 
than  1 /h . 


Proof  First  note  that 


,,,  N  -1  h 

(XVa  >  “  N 


N+Xh 


t — H_ 


*  •  • 


_ N_ 

N+Xh 


and  hence 


0  (XI  -a V1  « 
n  n 


(4.5) 


Now  let  z  e  Rn*f,+1)  be  given.  Then,  by  step  2,  the  function 


(t)  »  [ea  t  a  in  ] 


e  *nN  ,  t  >  o. 


is  square  integrable  on  the  interval  (0,»)  and  its  Fourier  transform 


*(iu)  «  ~~~  /  e  iU>t  z(t)  dt  •  — —  [(iwl  -aN)  1  *  I  ]  z(0) 


/2*  0 


/2* 


satisfies  lz> 
step  2  that 


,  u  =*  Izl  ,  .  Hence  if  follows  from  equation  (4.5)  and 

tr  L2(0,-,RnN) 


/  *  ■  N  /  II‘i»IM-aN)'1*  I J  ;<0„^  du 


nN 


2wh 


/  |z(ii»)  |  dm 

N  _ 
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I 


! 


2sh 

N 


J  l*(t)| 
0 


2 


nN 


dt 


2xh 

N 


2  -2£0t 

V  dt 


U(Q)I 


nN 


This  proves  step  3. 


Step  4  There  exists  a  constant  c  >  0  such  that  the  following  inequality  holds  for  every 
t  e  Rn(N+1>  and  every  N  > 

m 

j  l(iu>I-AN)_1zl2  du  4  2xc2IiI2. 

N  N 


Proof  Recall  that 


(1UI-AN)'1 


N  N  -1  N  T  N  N 

EiuiA  (ita,)  "V  Q  P 


+  T 


iw' 


u  e  r. 


(proposition  3.5).  By  step  3,  it  remains  to  establish  the  desired  inequality  for  the 

first  term  on  the  right  hand  side  of  this  equation.  Moreover,  it  follows  from  theorem  4.4 

that  the  operators  PN  are  uniformly  bounded  and  it  is  easy  to  see  that  the  operators 
M  M  •  N  T  N 

E^  and  (Ey)  *  (E Q  are  uniformly  bounded  on  the  imaginary  axis.  Thus  it  remains 

N  — 1 

to  prove  the  desired  estimate  for  the  term  A  (iu>)  .  But  for  |u>|  >  VAR  (n)  it  follows 

from  lemma  4.2  (ii)  that 

|4N(i<,r1|  -  I  l  (i*)-*'1  LN(i*>k1  4 

k*0 
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i 


This  inequality,  together  with  lemma  4.2  shows  that 

sup  /  |4N(i«)  1|2du  <  m. 

N>Nfl  - 

This  proves  step  4. 

Step  5  For  every  z  e  Rn^f*+'  *  and  every  N  > 

m  -j 

t  ■  A  .2  ,  2,  ,2 

J  le  zl  dt  <  c  III  < 

0  N  N 

Proof  Let  z  e  *n^N+1  '  and  z(t)  “  e*  fcz  for  t  >  0  and  N  1  NQ.  Then  z(t)  is 

square  integrable  on  (0,“)  and  its  Fourier  transform  is  given  by 
A  —  i  /  2  N  —  1 

z(lw)  *  <2r)  (iu)I-a  )  .  By  the  Fourier-Plancherel  theorem  and  step  4,  we  obtain 

m  N  oo 

/  le*  6  zl2  dt  -  (2W)'1  /  l(i#iI-AN)"1z|2  dm  <  c2|*l2  . 

o  -» 

This  proves  step  5  and  the  statement  of  the  theorem.  ■ 

REMARK  4.10 

The  uniform  exponential  decay  rate  -e  for  the  approximating  systems  (2.10),  (2.14)  which 

has  been  found  in  the  proof  of  theorem  4.9  is  always  larger  than  -1/h.  The  question 

remains  open  if  one  can  find  a  uniform  exponential  bound  for  the  approximating  systema 

with  the  exponential  decay  rate  »  t  c  where  <u  “  sup  {Re  HdetA(A)  “  0}  and  c  >  0 

o  u  jj 

N  A  t  N 

can  be  chosen  arbitrarily  small.  It  is  also  an  open  problem  if  the  operators  t  e  » 
converge  to  the  (compact)  operators  S(t)  e  L( M2 I  in  the  uniform  operator  topology  if 
t  >  h.  If  this  could  be  shown,  then  the  solution  to  the  uniform  stability  problem 
mentioned  above  would  be  an  immediate  consequence. 
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5.  CONCLUSIONS 


The  present  paper  studies  in  detail  certain  finite  dimensional  approximations  for 
linear  retarded  systems,  namely  the  averaging  approximation  scheme,  both  a  continuous  - 
and  a  discrete-time  version  as  well  as  the  relation  between  these  two*  It  turns  out  that 
these  finite  dimensional  approximations  show  -  under  several  aspects  -  precisely  the  same 
structure  as  the  underlying  RFDE.  In  particular,  the  duality  relations  are  of  the  same 
type  and  there  are  certain  structural  operators  which  play  an  important  role  for  the 
description  of  the  approximating  systems  and  are  analogous  to  those  which  have  recently 
been  introduced  by  BERNIER-MANITIUS  [3],  MANITIUS  [14),  DELFOUR-WANITIUS  [71  for  the  study 
of  RFDEs.  Moreover,  it  is  shown  that  these  operators  actually  converge  to  the 
corresponding  operators  in  the  theory  of  retarded  systems.  One  of  these  convergence 
results,  namely  theorem  4.6,  is  only  a  slight  extension  of  a  corresponding  result  by  REBER 
(17,  theorem  7.5). 

Based  on  this  detailed  analysis  of  the  structure  of  the  approximating  systems,  it  is 
shorn  that  both  the  discrete-  and  the  continuous-time  approximations  are  stable  in  a 
uniform  sense  if  the  underlying  RFDE  is  asymptotically  Btable.  Such  a  result  is  by  no 
means  obvious  and  not  all  approximation  schemes  have  this  property.  For  example,  it  is 
shown  in  KAPPEL-SALAMON  [12)  that  spline  approximations  for  RFDEb  can  never  have  the 
property  of  uniform  stability.  Nevertheless,  the  uniform  stability  result  provides  a 
crucial  step  in  the  convergence  proof  of  GIBSON  [9]  for  the  solutions  of  the  algebraic 
Rlccatl  equation.  Moreover,  the  structural  matrix  F  Introduced  in  this  paper  allows  a 
factorization  of  the  approximate  Rlccatl  operator  in  precisely  the  same  manner  as  it  is 
shown  in  KAPPEL-SALAMON  (12)  for  the  spline  approximation  scheme.  Finally,  it  seems 
likely  that  the  uniform  stability  results  of  this  paper  have  some  implications  for  the 
construction  of  finite  dimensl  nal  compensators  for  RFDEs.  This  is  a  research  problem  for 
future  investigations. 
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6.  APPENDIX 


In  this  section  ve  formulate  and  prove  two  general  functional  analytic  result  which 
are  needed  frequently  in  section  4. 


LEMMA  6. 1  Let  A  be  an  arbitrary  set  and  let  X,Y,z  be  Banach  spaces.  Moreover ,  let 
K°  e  L(X ,  Y ) ,  T°  e  L(Y,z),  T*  e  l(Y,Z),  a  e  A,  k  e  M,  be  bounded,  linear  operators  with 
the  propertlee 

(I)  ct ( K°x  |a  e  A,  x  e  X,  Ixl  <  1}  c  Y  is  compact 

( II )  T0y  «  lim  T^y  for  all  y  e  Y  uniformly  in  a  e  A. 

k*«  k 

a  a  a  a 

Then  T^K  tends  to  T  K  6  t(X,Z)  in  the  uniform  operator  topology  as  k  tends  to 
infinity  and  this  convergence  la  uniform  in  a  e  A. 


PROOF  Note  that  for  every  e  >  0  there  exist  finitely  many  x.,,...,  ex  and 

a ,,•••, a  e  A  such  that  for  every  a  e  A  and  every  x  e  X  with  Ixl  <  1  there  is  a 

n  a  01 4 

j  e  such  that  IK  x  -  X  *  e'  Hence  the  desired  uniform  convergence 

result  follows  from  the  Inequality 


IT®k“x  -  t'V’xI  <  [lT®l  +  lTai]lKax  -  X  Jx.l  ♦  It“k  jx  -  T°  K  Jx  1  . 
k  1  k  1  j  k  j  j 


The  next  result  is  a  quantitative  estimate  for  the  equivalence  of  L^-stabllity  and 
exponential  stability  for  strongly  continuous  semigroups.  This  equivalence  has  been 
proved  -  for  the  case  p  -  2  -  by  several  authors,  see  for  example  DATKO  [5),  CURTAIN- 
PRITCHARD  (41,  PRZYLUSKI  (161.  But  none  of  these  give  the  desired  quantitative  estimate 
which  is  essential  for  the  proof  of  uniform  stability  in  theorem  4.9.  Again  in  the  case 
p  »  2  such  a  quantitative  estimate  can  be  found  without  proof  in  GIBSON  (8).  We  mention 
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that  some  of  the  ideas  in  the  proof  of  the  theorem  below  are  taken  from  PRZYLUSKI  [8 


proposition  9)  and  ZABCZYK  [21,  theorem  5.1). 


THEOREM  6, 2  Let  S(t),  t  >  0,  be  a  strongly  continuous  semigroup  of  bounded 
operators  on  a  Banach  space  X  satisfying  the  exponential  bound 


IS(t)l 


L(X) 


<  Me 


u>t 


t  >  0, 


for  some  constants  M  >  1,  ui  >  0.  Moreover,  let  1  5  p  <  •  and  suppose  that 
exists  a  constant  c  >  0  such  that 


/  IS(t)xlp  dt  <  cP«xlP,  x  e  X. 
0 


Then,  for  every 


a  > 


P„P 
pc  M 


1 


there  exists  a  Y  “  Y(a,w,M,c,p)  >  1  such  that 


»S(t)ll(x)  <  Yeat,  t  >  0. 


REMARK  6.3  If 


pc 


J _ 

pmP 


<  a  <  0, 
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t 


linear 


(6.1) 


there 


(6.2) 


(6.3) 


(6.4) 


than  thara  exists  a  unique  T  >  0  satisfying 


e 


-apT 


+ 


T 


or  equivalently 


a 


cP^e^ 

T  ♦  cPm.^T  ' 


The  proof  of  theorem  6.2  shows  that  in  this  case  y 


y(a,w,M,c,p)  >  1  can  ce 


y 


Il-L 


cP^g^/P 

TVPC 


PROOF  OF  THEOREM  6.2 

Let  T  >  0  be  given  and  let  us  define 


€(T) 


T  ♦  cPMPe“PT 


0. 


Then  it  follows  from  (6.1)  and  (6.2)  that 


V  IS(T)kxlP 
k-0 


-  I  ~  I  IS( kT ) xlP  dt 
k-0  0 

-  T 

*  lxlP  +  l  -  /  IS(T-t)lP  IS(  (k-1  )T+t)xlP  dt 
k-1  T  0 


<  ,x,p  ♦  Z  [  sup 
(0 ,  T] 


IS(t)lP]  /  IS (t ) xlP  dt 
0 


lxlP 


(6.5) 


chosen  as 


(6.6) 


(6.7) 
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£_1lxlp 


for  every  x  e  X.  This  implies  that 


l  IS(T)k+1xlP  -  l  IS(T)kxlP  -  lxlP  <  <  1  — E )  l  IS(T)kxlP 
k-0  k-0  k-0 


and  hence 


IS(T)“xl  <  [  l  IS(T)k+mxlP]  Vp 


k-0 


<  <1-£)m/p[  l  IS(T)kxlP]Vp 
k-0 


<  (1-£)m/pe-1/plxl 


for  every  x  e  X  and  every  m  e  Now  let  t  -  nff  +  T  >  0  and  m  ( 

0  <  T  <  T.  Then  we  conclude  that 

IS( t) xl  <  IS(T)IIS(T)mxl 


<  M  e"  £_1/P  <1-t)B/plxl 


E  log(l-c)  iSSi^Sl  ,1.  ,, 
<  M  eUWt_1/PeP  p  T 


Ixl 


lpq( 1 -E ) 

lift  -1/p  -1/p  pT  !fflT  T 

Me  £  p  (1-c)  e  v  Ixl 


at.  . 

1  e  Ixl 
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M  and 


where 


a  -  a(T) 


1  ,  cP«Pe“pT 

• —  log  "•  -  — 

pT  T-tcpHpett*,T 


<  0 


and 


T  *c 

(compare  (6.S)  and  (6.6)).  Thua  the  statement  of  the  theorem  follows  from  the  fact  that 
a(T)  is  strictly  increasing  for  T  >  0  and  satisfies 


lim  a(T)  - - —  . 

T>0  pcPMP 
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